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Abstract 

Strong commensurate antiferromagnetism proximate to superconductivity is found in some members of the k-(ET) 
family, while a spin gap (SG) is found in the 6-(ET). Both k- and S-(ET) materials have frustrated triangular 
lattice structures. We show from calculations of spin-spin correlations within the effective half-filled band triangular 
lattice proposed for the k-ET, as well as for the real lattice, that long range AFM order is not obtained as a 
consequence of this frustration. We argue that some other mechanism reduces the magnetic frustration in these 
systems. We show that the low temperature magnetic states in these materials can only be understood if the 
effects of the cooperative charge and bond ordering transitions occurring at higher temperatures in these systems 
are taken into account. In the k-ET, this co-operative transition leads to unequal hole populations on the ET 
dimers that form the triangular lattice. 
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1. Introduction 

The family of BEDT-TTF (ET) salts k-(ET)2X 
feature the highest known Tc of the organic super- 
conductors. In K-(ET), antiferromagnetism (AFM) 
is found adjacent to superconductivity (SC) as in 
the high-Tc cuprate superconductors, leading to 
speculation that the occurrence of AFM and SC 
must be related to each other [1,2]. Within ex- 
isting theories of AFM and SC in the K-ET2X, 
these materials lie on both sides of a Mott-Hubbard 
metal-insulator transition, as described within a 
triangular lattice half-filled Hubbard Hamiltonian 
[1,3]. The schematic crystal structure of the k- 
ET lattice is shown in Fig. 1(a), where pairs of 
ET molecules with nearly parallel molecular planes 
form closely spaced dimers, with nearest neighbor 
dimers having nearly orthogonal orientations, and 
an average of one half hole per ET molecule. Be- 
cause the intradimer hopping integrals are much 
larger than the interdimer hopping integrals [4], 
the K lattice is an effective anisotropic triangu- 
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Fig. 1. (a) Schematic lattice structure of k-(ET), includ- 
ing the intermolecular hopping integrals bl, b2 and q [4]. 
(b) Anisotropic triangular dimer lattice; each dimer of (a) 
constitutes one lattice point in this effective p = 1 lattice. 
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lar lattice of Fig. 1(b), each site within which is 
one dimer of Fig. 1(a). The ratio t2/ti (note that 
here we denote the smaller frustrating bond as t2) 
is ^2/^1 = tii2l{tp + tq). Literature survey gives 
t^/ti = 0.5 - 0.7 for X=Cu[N(CN2]Br, 0.73 for 
X=Cu[N(CN2]Cl, 1.0- 1.10 for X=Cu2(CN)3, 0.55 
- 0.6 for X=l3 and 0.6 - 0.9 for X=Cu(SCN)2. The 
large ^2/^1 imply that frustration plays a signifi- 
cant role in these systems. 



Despite the anticipated large frustration, two 
important experimental features must be noted: 
(a) when AFM order does occur, it is commen- 
surate and (b) the (staggered) magnetic moment 
can be large. For example, the magnetic moment 
in K-(ET)2Cu[N(CN)2]Cl is 0.45 Ms/dimer [5]. 
This value is comparable to the magnetization 
per site in the isotropic square lattice Hcisen- 
berg spin Hamiltonian. Considering the triangular 
dimer lattice structure of the k-ET lattice, the 
experimentally observed magnetic moments arc 
inexplicably large, and we argue below that they 
cannot be understood within existing theoretical 
models. Mean field (Hartrec/Hartrec-Fock) calcu- 
lations are often cited as correctly reproducing the 
experimental magnitude of the AFM moment [3] . 
It is, however, well known that mean field theory 
greatly exaggerates broken symmetry in corre- 
lated systems: for example Hartree-Fock incor- 
rectly predicts long-range AFM in one dimension 
(ID), and gives qualitatively incorrect predictions 
for charge order (CO) in ID as well as 2D [6,7]. 

The anisotropic triangular lattice has been stud- 
ied by many authors, both within the Heisenberg 
spin Hamiltonian [8,9], and the half-filled Hubbard 
Hamiltonian [3,10]. Within the Heisenberg Hamil- 
tonian the antiferromagnetic phase can correspond 
to either the three-sublattice spiral phase or the 
two-sublatticc coUinear phase [1] . The spiral phase 
is commensurate only at a = J2/J1 = 1, where 
Ji ~ /f7 and J2 t%/U, respectively, and is pre- 
cluded from experiment. The calculated magnetic 
moment in the coUinear commensurate phase de- 
creases rapidly as a increases from zero, and van- 
ishes at a ^ 0.5 [8]. Estimates of a in k-(ET) 
range from 0.3 to 1.0. Only in the region of very 
small a are the calculated magnetic moments any- 
where near the experimental values of 0.3 0.4 
jiB- ^ recent Hubbard model calculation [10] on 
the anisotropic triangular lattice (Fig. 1(b)) using 
a path integral renormalization group approach for 
2 < U/ti < 10 found a nonmagnetic insulating 
state over a broad region of t2/t\. For example, 
at U/ti = 6, AFM is only present for t2/ti < 0.5 
(corresponding to a < 0.25), smaller than the es- 
timated a values for any k-(ET)2X [10]. 

From the above theoretical studies of the 
anisotropic triangular lattice it becomes clear that 
the effective half-filled model is not able to account 
for the large AFM moment observed in k-(ET). In 
this paper we investigate theoretically the ground 
state of the anisotropic triangular lattice, both in 
the effective half- filled limit (Fig. 1(b)). as well as 
the quarter- filled lattice (Fig. 1(a)). Our results 
confirm the previous result [10] that the k-(ET) 
system, considered from either the viewpoint of 




Fig. 2. Sublattice magnetization for the anisotropic trian- 
gular half-filled Hubbard model of Fig. 1(b) for 48 and 80 

site lattices. 

the effective half-filled band or the true quarter- 
filled model does not support an AFM state with 
large moments. We then ask the question as to 
whether the strong moment can be explained by 
the presence of an additional broken symmetry in 
the lattice of Fig. 1(a). Specifically, we consider 
the possibility that when creating one effective 
site from the two molecules in a dimer, the charge 
density on the two ET molecules within a dimer 
are not equal, i.e. there exists a CO state similar to 
that found in 0-{ET) [7]. We show below explicit 
calculations where the presence of CO does indeed 
increase the strength of the AFM. 



2. Results and Discussion 

2.1. Anisotropic triangular half- filled band 

We first present results for the half- filled effective 
model of Fig. 1(b). The Hamiltonian is 

- *2 X! (4,<TCi,'7 + cl.crCi.a + U ^ni,tn»,4, (1) 

where {ij) {[ij\) indicate nearest-neighbor bonds 
along the solid (dotted) lines of Fig. 1(b). The 
results shown were obtained for the ground state 
using the Constrained Path Monte Carlo (CPMC) 
method [11]. The CPMC method was checked 
against exact results for a small cluster (N=12 
sites) and found to be quite accurate. To quantify 
the strength of the AFM, we calculate the static 
spin structure factor, S{q). S{q) is peaked at g = 
(tt, tt) in the AFM state. The sublattice magnetiza- 
tion M may then be defined as = 3S'(7r, 7r)/7V. 
In the square-lattice Heisenberg model (U 00), 
M w 0.308 [12]. In Fig. 2 we plot M as a function 
of t2/t\. Morita et al found a M of ^ 0.1 at t2/t\ = 



0.5 and U = 6, and a transition to a nonmagnetic 

state at t^/ti ^ 0.6 [10]. Note that for these values 
of ^2/^1, Hartree-Fock erroneously predicts nearly 
saturated AFM moments [3]. While we have not 
yet performed detailed finite-size scaling of the 
CPMC calculations, our present results are con- 
sistent with reference [10], as the points plotted in 
Fig. 2 are upper bounds to the scaled values of M. 
K-(ET)2C1 with t2/ti ~ 0.73[4] clearly should be 
nonmagnetic within this model. 

2.2. Quarter-filled model: CO and AFM 

We now return to the quarter-filled lattice struc- 
ture shown in Fig. 1(a), within the following ex- 
tended Hubbard Hamiltonian: 

(ij) i,a 

We compare the spin-spin correlations of Eq. 1, 
with and without the presence of CO. Our goal is to 
establish that long range commensurate AFM can 
indeed occur in the K-lattice in the presence of CO. 

We first show exact diagonalization data for 
spin-spin correlations on the small quarter-filled 
lattice (16 sites and 8 holes) shown in Fig. 3. This 
lattice is equivalent to that of Fig. 1(a) with dimer 
bonds indicated by thick lines. Dimer 1 in Fig. 3 is 
connected to dimer 2 and 4 via effective ti bonds, 
and dimer 6 via effective t2 bonds. Superimposed 
on the lattice of Fig. 3 is a horizontal stripe CO 
of the type found in 9-(ET), adapted for the k- 
(ET) lattice. This CO state is parametrized by 
the site energy terms Ci^a- in Eq. 2. Note that in 
the actual systems CO is driven by V as well as 
electron-molecular vibration coupling, and the site 
energies in Eq. 2 merely model this mechanism. 
The coupled bond-charge-spin nature of the CO 
can then in principle lead to effective exchange in- 
tegrals between the dimers that are very different 
from those for the uniform system without CO. 
Prom charge considerations alone, once the CO 
insulating state is reached from the metallic state, 
there occur three different kinds of nearest neigh- 
bor intermolecular bonds, which we label as 1-1, 
1-0 and 0-0, where the numbers and 1 denote 
populations 0.5+5 and 0.5-6, respectively. The 
bond orders (c|^Cj,o- -|- h.c.) corresponding to the 
different kinds of bonds are different, and conse- 
quently the effective interdimer exchange integrals 
between the dimers of Fig. 1(a) in K-ET2X with 
CO are also different. The effective Heisenberg J 
between dimer 1 and dimer 2 is increased due to 
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Fig. 3. Horizontal stripe CO within the k-(ET) lattice, 
showing the spin ordering found in a 16-site exact calcu- 
lation. Bonds drawn with thick lines are dimer bonds; tq 
bonds are included in the calculation, but not drawn for 
clarity. The charge densities of the sites denoted by the 
unfilled circles are smaller than the ones denoted by the 
black circles in the CO state. Spin-spin correlations be- 
tween non-equivalent pairs of dimers are shown in the ta- 
ble, as a function of the CO strength, An = fi[„g — rismall, 
where ni,ig and rismall the large and small charge den- 
sities, respectively. Hamiltonian parameters are as in Fig. 
1(a), with U = 0.7 eV and V = 0.2 eV. AFM correlations 
are strengthened in the CO state. 

the stronger 1 1 bond, while the coupling between 
dimer 1 and dimer 4 is weakened, reducing the 
amount of frustration present, and strengthening 
the AFM ordering indicated by arrows. 

The table in Fig. 3 shows the exact spin-spin cor- 
relations, comparing uniform (e=0 in Eq. 2) and 
CO states. To enable comparison with the half- 
filled band, the correlations shown are spin-spin 
correlations between dimers, i.e. for dimers with 
sites (1,2) and (3,4), the dimer-dimer correlation is 
simply {SiS3) + {S1S4) + {S2S3) + {S2S4). Once 
CO is present, the spin-spin correlations clearly re- 
flect the AFM state indicated by the arrows, and 
increasing the CO strength continuously strength- 
ens the AFM. Note that the molecular sites with 
small charge density ( "0" above) do have spin den- 
sity, which is aligned ferromagnetically with the 
other site in the dimer. The weakness of antiferro- 
magnetic correlations for the uniform lattice with 
no CO is obvious, in agreement with the 1/2-filled 
calculation of Fig. 2. In the CO state, spin-spin 
correlations corresponding to the 1-0 intradimer 
bonds are necessarily ferromagnetic now (the same 
electron occupies both sites), but so are the spin- 
spin correlations corresponding to interdimei 1- 
bonds, thereby minimizing the spin frustration 
that characterizes the effective lattice of Fig. 1(b). 
Because of the large finite size effects associated 




Fig. 4. (a) Sublattice magnetization for quarter-filled lat- 
tices of Fig. 1(a), with (7 = 0.7 eV and V = eV. Hori- 
zontal axis indicates the amount of CO as in Fig. 3. Filled 
(open) circles indicate exact (CPMC) calculation for N = 8 
dimers. (b) Size dependence of for the quarter-filled 
lattice with no CO. The large-lattice limit of M vanishes. 

with our very small lattice, wc now turn to quan- 
tum Monte Carlo simulations of larger clusters. 



tained in either the quarter-filled or half-filled ef- 
fective models is much weaker than experimental 
values, and as frustration ^2/^1 increases from zero, 
the magnetization drops rapidly. Second, we have 
presented calculations assuming the presence of 
CO, and find that the magnetization increases con- 
tinuously as the CO is strengthened. The strength- 
ening of the AFM state occtirs due to reduced 
frustration in the presence of CO. Our result that 
strong AFM is absent in the K-lattice is in agree- 
ment with earlier work. One intriguing question 
that arises now is that experiments to date have 
not indicated CO in k salts. This could be because 
of our ignoring electron-lattice interactions in the 
Hamiltonian (2). As shown elsewhere, such inter- 
actions cause strong modulations of the hopping 
integrals in the presence of CO in the 1 /4-fillcd lat- 
tices [6,7]. It is possible that in the real materials, 
the diff'erences in CO are too small for easy detec- 
tion, but the resultant hopping integral modula- 
tions are large and reduce frustration. We are cur- 
rently pursuing self-consistent calculations of hop- 
ping integrals within model Hamiltonians that in- 
corporate electron-lattice interactions. 



2.3. Quarter-filled model: large-lattice results 

We have performed direct calculations on larger 
quarter-filled lattices (up to 160 lattice sites/80 
dimers), again using the CPMC method. For in- 
creased accuracy with the CPMC method, we have 
performed these calculations for F = 0. Results 
for these calculations arc shown in Fig. 4. Fig. 4(a) 
shows M vs. An for several lattice sizes (includ- 
ing a comparison with exact results for the 8 dimer 
lattice of Fig. 3). As in Fig. 3, we find that CO 
strengthens the AFM slightly. The = 8 results 
also indicate that CPMC may underestimate the 
strengthening of the AFM with charge order, in 
contrast to mean field calculations which overesti- 
mate the AFM. 

Fig. 4(b) shows the finite-size scaling of M for 
the uniform (no CO) quarter-filled lattice. We find 
that the scaled value of M is essentially zero, con- 
sistent with Fig. 2 and the results of Morita et al. 
[10]. 



3. Conclusions 

In conclusion, we have investigated the AFM 
state in the triangular lattice appropriate for k- 
(ET). We find two principle results: First, within 
the best estimate of model parameters appropriate 
for k;-(ET), the staggered magnetic moment ob- 
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